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ABSTRACT 

Aims. We study the effect of nonperturbative corrections associated with the behavior of particles after shell crossing 
on the matter power spectrum. We compare their amplitude with the perturbative terms that can be obtained within 
the fluid description of the system, to estimate the range of scales where such perturbative approaches are relevant. 
Methods. We use the simple Zeldovich dynamics as a benchmark, as it allows the exact computation of the full nonlinear 
power spectrum and of perturbative terms at all orders. Then, we introduce a "sticky model" that coincides with the 
Zeldovich dynamics before shell crossing but shows a different behavior afterwards. Thus, their power spectra only 
differ in their nonperturbative terms. We consider both the real-space and redshift-space power spectra. 
Results. We find that the potential of perturbative schemes is greater at higher redshift for a ACDM cosmology. For 
the real-space power spectrum, one can go up to order 66 of perturbation theory at z = 3, and to order 9 at z = 0, 
before the nonperturbative correction surpasses the perturbative correction of that order. This allows us to increase 
the upper bound on k where systematic theoretical predictions may be obtained by perturbative schemes, beyond the 
linear regime, by a factor ~ 26 at z = 3 and ~ 6.5 at z = 0. This provides a strong motivation to study perturbative 
resummation schemes, especially at high redshifts z > 1. 

In the context of cosmological reconstruction methods, the Monge-Ampere-Kantorovich scheme appears to be close to 
optimal at z = 0. There also seems to be little room for improvement over current reconstruction methods of the baryon 
acoustic oscillations at z = 0. This can be understood from the small number of perturbative terms that are relevant 
at z = 0, before nonperturbative corrections dominate. 

We also point out that the rise of the power spectrum on the transition scale to the nonlinear regime strongly depends 
on the behavior of the system after shell crossing. 

We find similar results for the redshift-space power spectrum, with characteristic wavenumbers that are shifted to lower 
values as redshift-space distortions amplify higher order terms of the perturbative expansions while decreasing the 
resummed nonlinear power at high k. 

Key words, gravitation; cosmology: theory - large-scale structure of Universe 
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The growth of large-scale structures in the Universe 
through the amplification of small primordial fluctuations 
by gravitat ional instabili ty is a key ingredient of modern 
cosmology (|Peeblesl 1 1 980t ) . and it can be used to constrain 
cosmological parameters through the dependence of the 
matter power spectrum on scale and redshift. On very large 
scales or at high redshifts, where the amplitude of the den- 
sity fluctuations is small, it is sufficient to use linear theory, 
whereas on small scales, in the highly nonlinear regime, one 
must use numerical simulati ons or phenomenologic al mod- 
els, such as the halo model (jCoorav fc Shethll2002h . which 
are also calibrated on simulations. In the weakly nonlin- 
ear regime one expects perturbative approaches to provide 
a useful tool, as they allow going beyond linear theory in 
a systematic and controlled fashion. S everal observational 
probes, such as we ak lensing surveys (jMassev et al.l 120071 
iMunshi et aDl2008h or measur e s of a coustic baryonic oscil- 
lations ( Eisenstein et al.l [l998, 2005), are mostly sensitive 
to these intermediate scales, and to meet the accuracy of 
future observations we need a theoretical accuracy of about 
1%. Phenomenological models typically have an accuracy 
of 10% in this range, while numerical simulations may suf- 



fer from finite resolution and finite size effects and require 
a long computational time if we need to obtain the power 
spectra over a fine grid of cosmological parameters. 

This has led to a renewed interest in perturbative ap- 
proaches, as it may be possible to improve over the standard 
pertu rbation theory (jGoroff et al.l Il986t iBernardeau et al.l 
2002) by using resummation schemes that allow system- 
atic partial resummat i ons of higher order terms . Thus, 
ICrocce fc Scoccimarrol (|2006bl lat) present a partial resum- 
mation of the diagrammatic series associated with the re- 
sponse function (propagator), within a high-fc limit, which 
provid es improved predictions for t he density power spec- 
trum (ICrocce fc Scoccimarrol 120081 ). On the other hand, 
IValageasT i 2007aD describes a path-integral formalism that 
allows applying the tools of field theory, such as large- TV ex- 
pansions, to compute the pow er spectrum an d higher order 
statistics like the bispectrum ((Valageas 2008) . One of these 
large- N expansions was recovered bv iTaruva fc Hiramatsul 
as a "closure theory" where one closes the hierar- 



chy of equations satisfied by the many body correlations 
at the third order, following the "direct i nteraction approx - 
imation" introduced in hydrodynamics ()Kraichnanl [1959) . 
This also improves the predictions for the power spectrum 
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on the scales probed by the baryonic acous t ic osc illations 
(jTaruva et all l2009t IValaeeas feNishimich.il I2010D. Other 
approaches have been proposed by Matarrese & Pictroni 
(l2007fl. using the dependence on a running high-fc cutoff, 
bv lPietronil (|2008l) . using a truncation o f the hierarchy sat - 
isfied by the many-body correlations, by Matsubara (2008), 
within a Lagrangian framework, and bv iMcDonaldl (|2007h . 
using a renormalization group technique. MostQ of these 
approaches start from the fluid description of the system, 
where the density and velocity fields obey hydrodynamical 
equations of motion. This corresponds to a single-stream 
approximation that neglects shell crossing. Then, the do- 
main of validity of most such perturbative schemes is lim- 
ited to wavenumbers where shell-crossing effects are negli- 
gible, even if we could sum all perturbative terms. 

This problem with the impact of shell crossing also ap- 
pears in the context of cosmological reconstruction, where 
one attempts to follow the matter distribution o bserved in 
a given galaxy survey back in time (Peebles 1989). Thus, an 
efficient algorithm for building such a reconstruction (which 
can then be used to estimate the velocity field) is provided 
by the Monge-Am pere-Kantorovich method, which neglect s 
multistreaming (jBrenier et al.ll2003tlMohavaee et al.l l2006'). 
Similar methods are used to reconstruct the baryon acous- 
tic oscillations (BAO) of the density power spectrum, in 
order to improve the accuracy of cosmological distance mea- 
surements and tighten the constraints on cosmological pa- 
rameters such as the amount and evolu tion of dark energy 
(jEisenstein et all 120071: ISeo et al.ll2010Fl . Then, by compar- 
ing the amplitude of nonperturbative corrections with the 
perturbative terms obtained within the fluid description, 
one can estimate the scale down to which these reconstruc- 
tion schemes can be used. 

Thus, to estimate the potential of such approaches, it is 
necessary to evaluate the effect of shell crossing on the mat- 
ter power spectrum. This question has already been inves- 
tigated by lAfshordil (|2007h by comparing the phenomeno- 
logical halo model with a modified variant where halos are 
collapsed to pointlike masses. In this paper we revisit this 
problem in a more systematic fashion, within the framework 
of the Zeldovich dynamics. Thus, we compare the Zeldovich 
dynamics with a second model (named the "sticky model" 
in the following), which only differs after shell crossing. 
Therefore, both power spectra have the same perturba- 
tive expansions and only differ by nonperturbative terms. 
Taking advantage of our being able to explicitly compute 
the full nonlinear spectra, as well as perturbative terms 
at all orders and these nonperturbative corrections, we can 
compare their respective amplitudes in detail. This allows a 
more detailed discussion of the importance of shell-crossing 
effects and of the scope of perturbative approaches. This 
also enables us to distinguish the dependence on the be- 
havior of the dynamics after shell crossing of the rise of the 
density power spectrum on mildly nonlinear scales, in the 
intermediate regime where the logarithmic power goes from 
A 2 (fc) ~ 1 to A 2 (fc) ~ 100. In addition, we can perform the 
same analysis for the redshift-space power spectrum, which 
is actually the quantity most directly observed in galaxy 
surveys. 



1 Tw o exceptions are the formalism developed in IVala gcas 
(2004), which applies to the Vlasov equation, and the 
Lagrangian approach of Matsubara (2008), which however is not 
valid beyond shell crossing. 



This article is organized as follows. We first recall in 
Sect. 12.11 the nonlinear real-space power spectrum asso- 
ciated with the Zeldovich dynamics, and its perturbative 
expansions in Sects. 12.21 and 12.31 We present our "sticky 
model" , which only differs from the Zeldovich dynamics af- 
ter shell crossing, in Sect. 12.41 and we give the associated 
nonperturbative correction. Then, we describe the numer- 
ical results obtained for a ACDM cosmology in Sect. [3J 
comparing the various perturbative and nonperturbative 
terms. Finally, we extend our analysis to the redshift-space 
power spectrum in Sect. 2J focusing on wavenumbers that 
are aligned with the line of sight. We conclude in Sect. [5j 



2. Computation of the density power spectrum in 
real space 

2.1. Nonlinear Zeldovich power spectrum 

As is well known, the Zeldovich approximation (IZeldovichl 
Il970f ) sets the Eulerian position, x(q, t), of the particle of 
Lagrangian coordinate q, equal to the position given by the 
linear displacement field, \tj,(q,t), 

x(q,i) = q + * L (q,t) with V q • *z = S L (q, t), (1) 

where Sl (q, t) is the linear growing mode of the density 
contrast, which is defined by 



P(x,t) 



P 



(2) 



Here p is the mean matter density of the Universe, 
we work in comoving coordinates and as usual we only 
consider the linear growing mode. It is well known 
(|Schneider fe Bartelmannl Il995t iTavlor fc Hamilton! fl996h 
that the explicit expression of the matter power spectrum 
can be derived from Eq.([T]) by using the conservation of 
matter, which reads as 



/?(x) dx = pdq whence 1 + <5(x) = 



dot 



()q 



(3) 



For an arbitrary displacement field VP, this also reads as 



l + «S(x) = J dq a D [x-q-¥(q)], 



(4) 



where 8d is the Dirac distribution, and this yields in Fourier 
space (for k ^ 0, that is, disregarding a term <5i)(k)): 



<f(k) 



dx 



,-ik : 



*(x) = 



dq _. 



k-(q+*) 



(5) 



(2tt) 3 v ' J (2tt) 3 

Defining the density power spectrum as 

(J(ki)S(k 2 )) = Mki + k 2 )P(fci), (6) 
we obtain from Eq.([5]), using statistical homogeneity, 



P{k) 



dq 

(2nf 



3 ik-[x(q)-x(0)]\ 



(7) 



Equation (J7J is quite general since we have not used the 
Zeldovich approximation ([1} yet. It shows how the density 
power spectrum is related to the statistical properties of 
the displacement field, for any mapping qi->x (which can 
include some shell-crossing, as in the Zeldovich case). The 
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great simplification provided by the Zeldovich approxima- Substituting this expansion (and the one associated with 
tion QJ is that in this case the quantity x(q) — x(0) is a the velocity field) into the equations of motion one obtains 



Gaussian random variable, so that the mean can be com 
puted at once as 



P(k) 



dq 

(27T)* 



e ik-q e -i((k-[* 1 (q)-* I .(0)]) 2 >^ 



Next, using the second relation ([T]) to compute the average 
((k • [^i(q) — S&l(0)]) 2 ), we obtain the explicit expression 



a recursion relation for the kernels F n , which allows terms of 
increasing order to be computed in a sequential manner. As 
is well known, for the Zeldovich dynamics, where \P = 
(g) we do need to follow this route, since by expanding the 
exponential |5) over and using *L(k) = i(k/fc 2 )^L(k), 
we obtain at once all terms, 



P(k) = 



dq 

(2tt)< 



e ikc V 



- J dw [1 — cos(w-q)] ■ 



F„(wi, .,w„ = — „— 

TJI Hit 



(18) 



■Pl(w) 



(9) 



It is convenient to perform the integration over the angles 
of the wav enumber w by expanding e lw q over spherical 
harmonics ([Schneider fc Bartelm ann 1995)), 

f dwe iw ' q QilJpl. P l (w) = k 2 I (q)+k 2 (l-3fi 2 )l2(q), (10) 
J w A 

where fi = (k • q)/(kq), and we introduced 
dwP L (w) je(qw), 



Then, substituting the expansion (1161) into the definition 
([6]) of the power spectrum and taking the Gaussian average, 
one obtains the standard perturbative series 



P{k) = J2 p(n) ( k ) with P {n) (k) oc (P L ) n . 



(19) 



n=l 



4tt 

««> - - 



(U) 



In particular, the two lowest order terms are 

P [1) {k) = P L {k), P {2) (k)=P 22 (k)+P 13 (k), (20) 

where P22 and Pi 3 arise from terms of the form (S^'S^) 
and (5^5 ( V), with 



where ji is the spherical Bessel function of order £. In par- 
ticular, the variance a 2 of the one-dimensional displacement ■ 22 

(k) 

field (or of the linear velocity field, up to a time-dependent 
multiplicative factor), reads as 



o(l*i! 2 > 



/o(0). 



(12) 



dwidw2 Sd(wi 

xP L ( Wl )P L (w 2 ), 
and 

P 13 (k) = -k 2 a 2 v P L (k). 



• W2 — k) 



(k-W 1 ) 2 (k-W 2 ) 



2w\w\ 



(21) 



(22) 



Therefore, the expression ([9]) also writes as 



P(k) = 



dq ,-fc 2 [ CT , 2 -/o(g)-(l-3 A1 2 )/ 2 (?)] 



cos(fcq/i) e 



In fact, for the Zeldovich dynamics it is not necessary to 
use this procedure, since by expanding the exponential in 
(13) the exact result ^ or (fT5|) we obtain all terms at once, 



f°° dq q 1 

Following [Schneider fc Bartelmann! (|l995h . we can perform P {n) {k) = j ^ ~ [ fc2 ( 7 o(<7) - 2/2(9) - OF 

the integration over the angles of q by expanding part of ~" ■■■ ' ' 

the exponential, and using the property 



dp, cos{kqp) (1 — p 2 ) 1 

which gives 

*dqq 2 u 



' kq 



U(kq), 



P(k) = 



k 2 K-I (q)+2I 2 (q)] 



2tt 2 



E 



Qkl 2 {q) 



(14) 

jt(kq)- 
(15) 



p=0 

/ 6fc/ 2 (g) ^ 
V 9 



jn- P [kq). 



(23) 



2.3. Renormalized perturbative expansion 

As pointed out bv lCrocce fc Scoccimarrol (|2006bh . it is use- 
ful to keep the q-independent exponential term e~ k CT " in 
Eq. ([T3")) and to define a "renormalized" perturbative expan- 



P(k) = er k2a ' P £\ k ) with P^(k) cx (P L ) 



(24) 

2.2. Standard perturbative expansion 

In the standard perturb ative approach (|Goroff et al.lll986t Indeed, while the standard perturbative terms P^ grow 
Bcr nardeau et all I2002T ) . one writes the nonlinear density increasingly fast at high k and show large cancellations, the 
contrast as a series over powers of the linear growing mode, "renormalized" terms e^^pj^ are positive and show a 
oo Gaussian decay at high k. Thus, each term peaks on a well- 

S(k) — S^(k) with 5^(k) oc (<5>l)™; (16) defined range of k and one can clearly see the contribution 

of diagrams of a given order to th e full nonlinear power 
spectru m. In fact, as noticed in ICrocce fc Scoccimarrol 
that is, ([2006bD . these terms also write as 



S^(k) = J dw 1 ..dw„fo(w 1 + .. + w„ -k)F„(wi,..,w„) Pj™}(fe) = n! J dw 1 ..dw„fc(w 1 + ..+w n -k) 
x S L (wx)..S L (-w n ). (17) x P n (wi,..,w„) 2 P L (w 1 )..P L (w;„). 



(25) 
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This can be seen at once by expanding the expression © 
and recognizing the square of the kernel F n given in Eq. (fTKl) . 

Thus, each new term pj™^ is associated with the kernel F n 
that couples n linear modes, so that one can follow the 
contribution of higher order mode couplings. Again, from 
Eq. ([T5l) we obtain all terms at once, 



A 2 n 1 

^£^ 2 («-2i 2 (*))F 

x (amy* jn _ p(kq)t (26) 

while the two lowest order terms are 

P%>(k)=P L (k), P<$(k) = P 22 (k), (27) 

where P 22 was given in Eq. (l21[) . 

For numerical purposes, it is convenient to obtain the 
standard perturbative terms (|19[) from the renormalized 
ones Expanding the prefactor in Eg. (|2~1T) gives 



71 — 1 -. 

pM(k) = y / -r(-k 2 alrP^(k). 



^ pi 



(28) 



On the other hand, to compute the full nonlinear power 
spectrum (TTSl) it is convenient to subtract the two terms 



e k ' a ^[P^J (fc) + Pa v '(k)], so as to avoid integrating over 
badly behaving terms for q — > 00. 



2.4. Nonperturbative correction 

The expressions recalled in the previous sections give the 
nonlinear power spectrum and its perturbative expan- 
sions for the usual Zeldovich dynamics, where particles fol- 
low the linear trajectories ([I}. This includes shell cross- 
ing and leads to a decay of the nonlinear power spec- 
trum at high fc, as particles freely escape to infinity, 
and these random trajectories erase small-scale features 
([Schneider fe Bartelmarml Il995t iTavlor fc Hamilton! fl99l 
I Valageas! l2007bt iBernardeau fe Valageadl2010air However, 
because of the simple nature of the Zeldovich dynamics, the 
full nonlinear result (fT3")) can be obtained by resummin^l 
the perturbative expansions (fT9| or ((24]) . The latter can be 
obtained from the usual hydrodynamical equations of mo- 
tion in Eulerian space, which actually break down at shell 
crossing. 

Our goal in this work is to compare the perturbative 
terms (flU|) and (|2~4"|) with the nonperturbative terms that 
are associated with the physics that takes place after shell 
crossing. Within the framework of the Zeldovich dynamics 
considered in this article, this means that we wish to com- 
pare the previous results with those that would be obtained 
for a second dynamics, which coincides with the Zeldovich 
dynamics until shell crossing. Then, the difference between 
both power spectra would give us the amplitude of the cor- 
rection due to shell-crossing effects, which are generically 
nonperturbative. 

To introduce this second model, let us first recall that 
the linear growing mode of the displacement field, \I>L(q), 



2 The radius of convergence of such series can be finite (e.g., 
for a power-law line ar spectrum Ph(k ) oc fc' 1 with n — —2) or 
zero (if n < -2), see I Valageas! (|2007bh . 



and the asso ciated peculia r velocity field, v^(q) oc ^(q), 
are curlfree (|PeeblesllT98 Cfi. and are derived from a veloc- 
ity potential, vj, = — V • \l- Using the Poisson equation 
we can see that this potential is equal to the linear gravita- 
tional potential up 1 to a time-dependent factor (jPeebles! 
1980; Vcrgass ola et al.lfl994l) . In particular, the Lagrangian- 
space to Eulerian-space mapping defined by the Zeldovich 
dynamics ([!]) derives from a Lagrangian potential cp£,(q), 



x(q,t) 



with 



fL{<\,t) 



dq 



a(t) 



= *i(q,*), 



(29) 



(30) 



where is the linear gravitational potential, which obeys 
the Poisson equation, 



a(t) 



(31) 



Here Q is Newton's constant and a(t) the scale factor. Then, 
from Eq.Q the nonlinear density contrast is given by the 
Hessian determinant of pl , 



1 + a(x) 



dot 



9x 
dq 



det 



\dq l dq J 



(32) 



At early times, t — > 0, (and if there is not too much power 
on small scales), the Lagrangian potential tpL is dominated 
by the first term in Ea. (l3Ul) and x(q, t) — > q. Then, the 
Hessian matrix (d 2 >pL/dqidqj) is definite positive and goes 
to the identity matrix, whence S — ¥ 0. As time increases 
and structures form, the Lagrangian potential becomes in- 
creasingly sensitive to the fluctuations in the linear grav- 
itational potential, and the Hessian determinant deviates 
from unity. However, the eigenvalues remain strictly posi- 
tive until shell crossing, which means that the Lagrangian 
potential ph remains a strictly convex function. Then, at 
shell crossing one eigenvalue goes through zero and becomes 
negative (generically the collapse proceeds at different rates 
along the three axes). Thus, as is known (IVergassola et ail 
119941: iBrenier et al.l [2?5ol IBernardeau fc Valageasll2010bl ). 
the onset of shell crossing is associated with the change in 
sign of the Hessian determinant of <pl and with the loss 
of convexity of the Lagrangian potential p^ (the Hessian 
matrix is no longer positive-definite). 

A well-studied dynamics that agrees with the Zeldovich 
dynamics un til shell crossing is p rovide d by the "adhe- 
sion model" (|Gurbatov et al.lll989L Il99lt IVergassola et al.l 
I1994D. More precisely, the "geometrical adhesion model" 
(IBernardeau fc Valageas! l2010bl : IValaeeas fc Bernardeaul 
I2010TI leads to replacing the linear Lagrangian potential 
p>L{q), which defines the mapping q i-> x through (|2^|) . 
by a nonlinear Lagrangian potential p given by the convex 
hull of pl, that is, p = conv(</?£). Then, particles no longer 
cross each other but form shocks. Moreover, this second dy- 
namics still coincides with the Zeldovich dynamics outside 
of shocks, where ph coincides with its convex hull. 

In this article, we do not compute the density power 
spectrum associated with this "adhesion model", or an- 
other explicit dynamical system, which is a difficult task. 
Since we are merely interested in the amplitude of the ef- 
fects associated with shell crossing, we simply make use of 
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the property that shell crossing is associated with the loss 
of convexity of the Lagrangian potential y>/,(q). Then, we 
note that as long as <pz,(q) is strictly convex its restric- 
tion along a line that goes through two arbitrary points 
q^ and q^ is also str i ctly c o nvex, which implies (see also 
iNoullez fe Vergassolal (|1994l ); iBrenier et all (|2003D L 

before shell crossing, for any qa ^ <Ib ■ 

[x(q s ) - x(cu)] • (ciB ~ <U) > 0. (33) 

Indeed, we can choose a coordinate system so that = 
and qs = (9,0,0), and from strict convexity we obtain 
d 2 (pL/dqf > 0, whence dx\/dq\ > along the first axis and 
xi (q) — xi(0) > 0, which reads as ([3"3"]) in its general fornfl 
This means that the projection of the Eulerian separation 
vector, Ax = xg - x^, onto the Lagrangian vector, Aq = 
qs q.4, is positive, so that Ax lies in the forward half- 
space delimited by the plane orthogonal to Aq. 

Going back to the general expression ([7]), the density 
power spectrum is fully determined by the mean (e lk Ax ) q , 
for any q, where we note Ax = x(q) — x(0). To compute 
this quantity we can take q along the first axis, that is, 
q = qei with q = |q|. Then, before shell crossing we have 
from the constraint (|3"3"|) the property Axi > 0. Therefore, 
we consider the "sticky model" defined by: 

"sticky model", for q = |q| ei : Axi — max(Axz,i, 0), 

Ax 2 = Ax L2 , Ax 3 = Ax L3 , (34) 

where Ax^ is the linear Eulerian separation, as given 
by Eq.([T]). Thus, this second model only differs from the 
Zeldovich dynamics when the parallel linear Eulerian sep- 
aration, Axli, is negative, in which case we set it equal 
to zero. This is thus a simplified version of the "adhesion 
model" , as once Axli reaches zero, it remains equal to 
zero forever. However, the model ([34]) cannot be explic- 
itly derived from the "adhesion model" , since we take nei- 
ther transverse directions nor larger scales into account. 
Therefore, we use the more generic name "sticky model" , 
to refer to this sticking along one direction for the pair sep- 
aration. 

It is clear that the condition Axli < 0, where the two 
models differ, is only a sufficient condition for shell cross- 
ing, but it is not necessary. Thus, it is a local condition 
that does not take the "cloud-in-cloud" problem into ac- 
count : even though no shell crossing seems to have ap- 
peared on scale q yet, it may happen that this region is 
enclosed within a larger domain of size L that has al- 
ready collapsed, so that particles in t he smaller domai n 
have actually experienced shell crossing (iBond et al.lll991l) . 
In terms of the Lagrangian potential ifL, which defines the 
Lagrangian mapping, qn-x, through Eq. (j2"9"l) . the absence 
of shell crossing on a small domain of size q means that lqj 
is equal to its convex hull in this domain (|Vergassola et all 



3 In ID the property (|33[) is an obvious consequence of the 
absence of shell crossing, but this is not the case in higher di- 
mensions. Indeed, the result (1331) makes use of the constraints as- 
sociated with the fact that the Lagrangian mapping derives from 
a Lagrangian potential as in (|29[) . This relies on the Zeldovich 
dynamics and on the linear growing mode of the velocity and dis- 
placement fields being curlfree. For the gravitational dynamics, 
there is no such clear signature of the absence of shell crossing, 
since even in regular regions the displacement field ^(q, t) de- 
velops rotational terms at the third order of p erturbation theory 
l|Bucherdll99l : iBernardeau fc Val ageas 200|). 



119941: iBec fc Khanin1l2r)07t IBernardeau fc ValageaslEoiObfl . 
However, the construction of the convex hull is a global 
problem, as one must consider the behavior of ifL (q) over 
all the space, thereby taking into account the "cloud-in- 
cloud" problem, while in the definition of the model (fMf 
we only check a weaker condition, since we only consider 
the two points and q. This means that we somewhat un- 
derestimate the effects of shell crossing, but we can expect 
to obtain a reasonable estimate of their amplitude because 
the probability of collapse decreases on larger scales and 
we perform a statistical integration over the angles of q in 
Eq.0. 

From the previous discussions, the "sticky model" ([Mf 
and the Zeldovich dynamics ([I} coincide before shell cross- 
ing, since then we have Axli > 0. This implies that both 
theories coincide at all orders of the perturbation theory; 
that is, they show the same expansions (|23l) and (j2"6"]l over 
powers of Pl- However, they differ through nonperturba- 
tive terms, which arise from their different behaviors after 
shell crossing. 

From Eq.(j33]) we obtain A^i = max(A* L1 , -q), for 
q = (J 1 ex, and 

^ik-Ax^ = e ikc t e ~I dw[l-cos(wq)] |k J Pl(w) 

xe i*M P d^e-tA^/K^A^ (35) 

J-oo V27TCT|| 

where we factorized the result associated with the usual 
Zeldovich dynamics © in the first two terms and we intro- 
duced the variance of the linear longitudinal displacement, 

of(q) = ((A* L1 ) 2 )=2 J dw[l-cos( Wl g)]^|p L H (36) 

= 2a 2 v -2I Q (q)+AI 2 (q). (37) 

Separating the contribution from A^^i < —q, we obtain 
for the last two terms of Eq. (|55|) , 



**;*if/*UtiA* lN = 1 + e ifeH / dA ^Ll p -(A^ 1 ) 3 /(2^) 



x ( e - ifc i9 - gifeiA*i 



-1 i 1 ifc?a"n— ifciQ r 

= 1+ - e 2 1 ii H erfc 



2ito\\ 



(38) 

_ , -^erfc = " ,(39) 



q 



where erfc(z) is the complementary error function (ex- 
tended to the complex plane), 



erfc(z) = I dt e"' 



(40) 



Then, substituting into Eas. (|3"5l) and J?]), we can see 
that the density power spectrum of the "sticky model", 
Psticky(fc), is equal to the usual Zeldovich power spectrum 
obtained in Sect. I2.1l plus a correction term P s . c .(fc), 

Pticky(fc) = Pzcl(fc) + Ps.c.(fc), (41) 

with 

P sc (k) = - f dq c -k 2 (i-p?)W?.-i n M-i2M] c -g 2 /(2^ 
2 J (27r) 3 



x < w 



i<7 



y/2a\\ 



iq — k\i(j\ 



(42) 
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Here we introduced Faddeeva's function, 
w(z) = e~ z erfc(— iz), 



(43) 



which satisfies the asymptotic expansion 

(jAbramowitz fc Stegunlll97fl 



37T 

| arg(z)| < — , z — > oo : w(iz) 



— ( - — 

FUz { ~ 2z 2 



.(44) 



Expression (|42|) also reads as 



Ps c (fc) = r e -o 2 /^D f 1 dM e -^-^)K-io-h\ 

( 27r ) Jo 



xRe < w 



V2m\ 



iq — k\xa\ 
V2a\ 



(45) 



As expected, we can check on Eas. ([4*2)1 and (|4"5)) . using the 
behavior (|4"4")) , that the correction due to shell-crossing ef- 
fects is nonperturbative, in the sense that, because of the 
, its expansion over powers of the amplitude 



term e 



-97(2^) 



of the linear power spectrum is identically zercQ. 

At low k the shell-crossing contribution (|45p behaves as 
-P s c .(fc) oc k 2 . This agrees with the generic beha vioiB as- 
socia ted with small-scale redistributions of matter (|Peebles! 
Il974h . Thus, the "sticky model" provides an explicit exam- 
ple of a nonperturbative power spectrum, which satisfies 
these generic behaviors. 

3. Numerical results in real space 

We now describe the numerical results we obtain for the 
Zeldovich and "sticky model" power spectra, Eqs. tfTB")) and 
(|41[) . as well as their common perturbative expansions. We 
consider a ACDM cosmology, with fi m = 0.279, ft\ = 
0.721, cr 8 = 0.817, and h = 0.701. 



3.1. Logarithmic power 

We show in Fig. [T]the power per logarithmic interval of k, 
defined as 



A 2 (fc) = 47rfc 3 P(fc), 



(46) 



for redshifts z = and z = 2. We compare the Zeldovich 
nonlinear power spectrum (I13|) with its perturbative ex- 
pansions (|T9"]) and (|2~4"]) and the nonperturbative correction 
(|45[) . As is well known, higher order terms of the standard 



4 This is not a large-scale or early-time expansion, which 
would require knowledge of the shape of Pt(fc). As for the stan- 
dard perturbative expansions of Sects. 12.21 and 12.31 one intro- 
duces for instance an auxiliary parameter A, such as 5l —¥ \5l 
or Pl —¥ X 2 Pl, and looks for an expansion over powers of A. 
Since a 2 oc q 2 for g — > if Pt(fe) decreases faster than k~ s at 
high k, as for CDM power spectra, the integration over q does 

not transform the exponential factor e 9 ^ 2<T ll- ) into anomalous 
power laws, and the expansion of P s . c . (k) over powers of the am- 
plitude of the linear power spectrum is indeed identically zero, 
hence nonperturbative. 

5 As described in Peebles (1974), considering momentum con- 
servation, in addition to matter conservation, leads to a k 4 tail 
at low k. However, since the "sticky model" (|34[) does not satisfy 
momentum conservation but only matter conservation, it gives 
rise to the k 2 tail. 



perturbative expansion (fl"9")) grow increasingly fast at high 
k with changes in sign and large cancellations between dif- 
ferent orders. The "renormalized" perturbative expansion 
AMI) gives posi tive terms (see Eq.(|2"5"|)) that p eak on a well- 
defined range (|Crocce fc Scoccimarroll2006bl) and are much 
easier to distinguish. Thus, while we only plot the first three 
orders of the standard expansion, P« = p L and the ab- 
solute values |P( 2 )| and l-P^I, we plot the first five orders 

of the "renormalized" expansion, 

pW t0 as we n 

as 

p£°>, P^\ pj 20 >, pjf, Pif) and PJI 0) (always multi- 
plied by the factor 47rfc 3 of Ea. (j4"6"[) '). As we go to higher 
orders, contributions become narrower and more densely 
packed, which implies that as we go deeper in the nonlin- 
ear regime, we need increasingly more perturbative terms 
per logarithmic interval of wavenumber. 

We also plot a few partial sums of expansion (f2"4"]l , that 



is, e~ k '< Yln=i p °'\ k )i with N = 5 > 10 > 15 > 2 °. and 30 - We 
can check that they agree with the full nonlinear power 
spectrum (|13[) until the wavenumber associated with the 

peak of e~ k ""Pi^, after which they follow the Gaussian 

,2 2 

decay associated with the prefactor e <J ". These partial 
sums are also slightly more efficient than those obtained 
from the standard expansion (|19p at the same order (not 
shown in the figure). 

As expected, the nonperturbative correction (1451) is very 
small on quasi-linear scales, so that there is indeed a range 
where higher order terms of the perturbative expansions 
([1"9"| and (|2"4"|) (i.e. beyond the first order associated with 
the linear regime) are relevant. At higher fc, the nonpertur- 
bative correction becomes dominant and the perturbative 
expansions become irrelevant, since one can no longer ne- 
glect the physics beyond shell crossing. 

In the highly nonlinear regime, we recover the 
well-known decay of th e Zeldovich logarithmic non- 
linear power spectrum ([Schneider & Bar telmannl 119951 : 
iTavlor fc Hamilton! Il996l IValageasll2007blk This isbecause 
particles escape to infinity after shell crossing (i.e. keep 
moving on their straight trajectories), and these random 
trajectories (in the sense of random initial conditions) erase 
small-scale features in the density field. This is expressed 
by the Gaussian decaying factor e~ k [ <T i.~ / o-( 1 - 3 A' K2] m 
Eq.(|13p. where the quantity within brackets is always pos- 
itive as can be seen from expression ([5]). Even though this 
leads to a decay at high k, the falloff is not Gaussian be- 
cause of the integration over the Lagrangian distance q 
(jTavlor &: Hamiltonlll996l ). In particular, for power- law ini- 
tial power spectra, Ph{k) oc k n with — 3 < n < —1, one 
finds that A 2 (fc) oc fc3(n+3)/(n+i) (jvalaeeasl [2007bh (the 
Zeldovich dynamics is not well-defined for n > — 1 because 
of ultraviolet divergences). Contrary to the gravitational 
case, the nonlinear power spectrum decreases faster at high 
k for higher values of n. This is due to the greater smearing 
out of small-scale features by the larger amplitude of the 
random linear displacements at small wavelengths. 

It is interesting to note that, thanks to its shell-crossing 
correction (|45[) . the nonlinear power spectrum (I41[) of the 
"sticky model" does not show this fast decay, and its loga- 
rithmic power still increases in the nonlinear regime. This 
is due to the prescription (|34l) . which in a sense prevents 
particles from escaping to infinity in one direction, as they 
stick together. This approximately models the formation of 
Zeldovich pancakes, associated with collapse along one axis 



P. Valageas: Impact of shell crossing and scope of perturbative approaches in real and redshift space. 



7 




Fig. 1. The power per logarithmic interval of k, as defined in Eq. (|4"oT) . at redshifts z = (left panel) and z = 2 (right 
panel). The solid lines are the linear power spectrum "Pj,", the nonlinear Zeldovich power spectrum "Pz e i", given by 
Eq. (fT5|) . the nonpcrturbative correction "P s . c .", given by Eq. (|4"5)) , and the sum Pz c \ + P s . c . associated with the "sticky 
model", Eq. (|4"Tj) . The dot-dashed lines, which grow very fast and are only partly drawn, are the absolute values of the 
standard perturbative terms \P^ \ and |P ( - 3 ''|, from Eq. (j2"3")) . The lower dashed lines are the "renormalized" perturbative 

terms e'^'pi"' of the expansion l[2l|), from Eq.®, for n = 1 to 5 and n = 10, 15, 20, 30, 50, and 70. The peak moves to 
higher k as the order n increases. The dotted lines that follow the nonlinear Zeldovich power spectrum until a Gaussian 

decay are the partial sums e~ k a » X^=i Parity, with TV = 5, 10, 15, 20, and 30. All terms are multiplied by the factor 
47rfc 3 of Eq.®. 



(although there is no precise relationship, since Eq. ([34]) is 
only a statistical model and does not consider the "cloud-in- 
cloud" problem). Thus, small-scale structures are no longer 
erased, since we keep a trace of planar features. This is ex- 
pressed by the factor (1 — /i 2 ) in the exponential argument 
in expression (|4"5)l . which suppresses the Gaussian decay- 
ing term of the form e for \i ~ 1. This gives a width 
Afi ~ k~ 2 , hence P s . c .(fc) ~ k~ 2 at high k, as would be the 
case for a density field where planar objects are the rele- 
vant nonlinear structures (i.e. bi-dimensional structures, as 
opposed to pointlike masses or lines for instance). Contrary 
to the Zeldovich power spectrum there is no dependence on 
n for the high-A; slope, for power-law linear power spectra 
Pt(fc) cx k n with —3 < n < —1. This gives the univer- 
sal asymptote Ag ticky (fc) cx k for the nonlinear logarithmic 
power of the "sticky model" at high k. 

3.2. Perturbative and nonperturbative contributions 

To clearly see the range of scales where perturbative 
schemes are relevant, we plot in Fig. [5] the ratios of the suc- 
cessive terms of the "renormalized" perturbative expansion 
|24|) with respect to the nonlinear power spectrum of the 

"sticky model", e'^^ P^\k)/ P sticky {k), for n = 1 to 5, as 
well as n = 10, 15, 20, 30, 50, and 70. We also plot the ratio 
-Ps.c. (k)/P 3 ticky(k), to compare with the amplitude of the 
nonperturbative correction associated with shell-crossing 
effects. We consider the four redshifts z = 0,1,2, and 3. 
Here we focus on the "renormalized" perturbative expan- 
sion (|24|). rather than on the standard expansion (fl9]) . to 
avoid the interferences brought by the changes in sign and 
the cancellations between various terms. 



We can see that the potential of perturbative expansions 
grows at higher redshift for a ACDM power spectrum, be- 
cause of the change in slope of Pl(&). Thus, at z = 3 we 
can go up to the order n = 66 before the perturbative 
term becomes subdominant with respect to the nonper- 
turbative correction, while at z = the crossover takes 
place at n = 9, as the term e~' £ " cr « P<r„ (&) is the lowest 
order one that is fully below the nonperturbative correc- 
tion P s .c. (fc). This agrees with the observation that in the 
gravitational case perturbative schemes (an d resummation 
appro aches) seem to fare better at higher z (ICarlson et al.l 
2009). We can see that perturbative schemes are relevant 
over roughly one decade over k. Thus, if we require an ac- 
curacy of 1% at z = 0, linear theory is sufficient up to 
fciioop ~ 0.033ft. Mpc" 1 , while higher order perturbative 
terms allow reaching fc s c . ~ 0.23ft Mpc -1 . At higher k one 
must take the nonperturbative correction associated with 
shell-crossing effects into account, which implies going be- 
yond the fluid approximation and requires new approaches. 

To help the reader, we give in Table [T] the wavenumber 
^lioop) below which the linear term is enough to reach a 
1% or 10% accuracy, for the four redshifts shown in Fig. O 
We also give the wavenumber k s . c . beyond which the non- 
perturbative correction is required to reach an accuracy of 
1%, 10%, or 50% (in units of P st icky(fc)). Thus, the inter- 
val [fciioopj k s . c ] gives the range where perturbation theo- 
ries based on the fluid description are relevant. Of course, 
this range shifts to higher k at higher redshift. It is in- 
teresting to note that this range is also broader at higher 
redshift, as the slope of the CDM linear power spectrum on 
the relevant scales changes slowly. The last column gives 
the last order, n s c ., of the "renormalized" perturbative ex- 
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Fig. 2. The dashed lines are the ratios of the "renormalized" perturbative terms of the expansion (|24|) to the nonlinear 
"sticky model" power spectrum, e~ k * a * P^\k) / P a tvdky{k) , for n = 1 to 5 and n = 10,15,20,30,50, and 70. The solid 
line is the ratio of the nonperturbative correction to the nonlinear power spectrum, -P s .c.(fc)/-Psticky(fc)- The four panels 
correspond to redshifts z = 0, 1, 2, and 3 (with a shift of the plot towards higher k at higher z). 



pansion that is not fully below the nonperturbative term. 
As noticed above in Fig. [5J this expansion order is sig- 
nificantly higher at higher redshift. This corresponds to a 
greater potential for perturbative schemes. However, n s , c , 
grows faster than the logarithmic width of the perturbative 
range, [lnfcn oop , lnfc s . c .]. Indeed, as shown in Fig. [3J peaks 
associated with higher order perturbative terms are increas- 
ingly narrow on the In k axis. This implies that to multiply 
the upper wavenumber k, defined by a fixed accuracy, by 
a given amount, one needs to add an increasingly greater 
number of new perturbative terms. 

In practice, we do not expect that perturbative terms 
will be computed up to such high orders, since in the case of 
the gravitational dynamics this would involve multidimen- 
sional integrals that are beyond the reach of current nu- 
merical possibilities. However, resummation schemes allow 
one to consider parts of such higher order terms (actually, 
an infinite number of diagrams that contribute to terms of 
all orders). Then, the hope is that such methods can effi- 
ciently sum most of the contributions of higher order terms 
and accelerate the convergence on weakly nonlinear scales. 
The comparison displayed in Fig. [5] shows the potential of 
such methods (i.e. the best result one can obtain, for the 
"sticky model" considered here), which appears to be quite 
significant at higher redshifts, z > 1. 



Since the Zeldovich dynamics provides a reasonable 
approximation of the gravitational dynamics down to 
weak ly nonlinear scales (jColes et al.lll993t iPauls fc Melotd 
1995), and its accuracy can be improved by imple- 
menting the "adhesion model" that only differs after 
shell crossing ([Weinberg fc Gunnlll990l iMelott et al.lll994t 
ISathvaprakash et al.lll995[ ). we can expect that to a large 
extent these results still apply to the gravitational case. 
In particular, while we find that at z — the nonpertur- 
bative correction to the density power spectrum is around 
1% at k ~ 0.2 3/1 Mpc" 1 and 50% at k ~ 0.9h Mpc" 1 , 
lAfshordil (|2007fl finds the wavenumbers k ~ O.lh Mpc -1 
and k ~ 0.85h Mpc -1 with a phenomenological "sticky 
halo model" . It is comforting that these two very different 
approaches give similar results. Then, the property that the 
range where perturbative schemes are relevant is greater at 
z = 3 than at z = 0, with a higher order n s c ., should remain 
valid. 

This is confirmed by the analysis of unified mod- 
els that combine perturbation theories with halo mod- 
els, for the gravitat i onal case. Thus, as shown in 
IValaeeas fc Nishimichil (|2010h . at higher redshift the in- 
termediate range, where the power spectrum departs from 
one-loop perturbation theory but is not yet well described 
by the "one-halo" contribution, becomes wider. There, the 
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1% 


0.033 


0.23 







10% 


0.082 


0.45 


9 




50% 




0.9 






1% 


0.043 


0.44 




1 


10% 


0.11 


1.1 


18 




50% 




2.2 






1% 


0.057 


1.2 




2 


10% 


0.14 


2.3 


37 




50% 




6.4 






1% 


0.07 


2.2 




3 


10% 


0.18 


5.2 


66 




50% 




10.4 





Table 1. The dependence on redshift z of the wavenum- 

bers where the ratios e~ fc2<T "Pi^(fc)/P st i C ky(fc) (at fcn oop ) 
and Ps.c.( k ) / Psticky(k) (at fc s . c .) reach 1%, 10%, and 50%. 
The last column gives the last order n s . c . of the "renor- 
malized" perturbative expansion that is not fully below the 
nonperturbative term. 



"one-halo" contribution plays the role of the nonperturba- 
tive contribution (1451) , as it is also fully nonperturbative and 
decays at low k as k 2 (if w e only consider matter conserv a- 
tion, see the discussion in IValaeeas fc Nishimichil (|2010r ib 
Then, the broadening of this intermediate range also sug- 
gests that higher orders of perturbation theory become rel- 
evant, as explicitly found in Figs. Q] and [5] 

In the context of cosmological r econstruction^!, a com - 
parison with iV-body simulations (Mohava ee et al.l [2006) 
shows that the Monge-Ampere-Kantorovich method is able 
to recover the nonlinear displacement field down to ~ 3/i _1 
Mpc at z = 0, which corresponds roughly to k ~ 2h Mpc -1 . 
We can see from Table [T] that this is a very good result, as 
going to smaller scales requires taking shell crossing into 
account (in fact, at k — 2h Mpc -1 nonperturbative correc- 
tions have already started to dominate). Thus, the Monge- 
Ampere-Kantorovich method appears to be close to opti- 
mal at z = 0, because it goes as far as any scheme that 
disregards shell-crossing effects can be expected to go. This 

6 There, the problem is to reconstruct the past dynamical 
history of a given region of the sky from the knowledge of its 
present density field. A key observation is that in the linear grow- 
ing mode, the velocity field is curlfree and related to the linear 
gravitational potential by a relation of the form v_l oc — \7<&l, 
while the linear density field is given by the Poisson equation, 
8l oc A^l- Then, the number of unknowns (the linear grav- 
itational potential) is equal to the number of constraints (the 
present nonlinear density field), so that it makes sense to look 
for a reconstruction defined by the nonlinear density field alone 
(but this does not ensure uniqueness). The Monge-Ampere- 
Kantorovich method makes use of an additional assumption, 
that the Lagrangian map, q \— > x, can be derived from a con- 
vex potential <p(q), as in Eq. (|29p . to derive a unique solution to 
the "displacement reconstruction", i.e. the inverse map x H y q 
ijBrenier et al.| [2003; Mohavaee ct al. 2006). Such methods can- 
not describe multi-streaming, and by estimating the scale where 
shell crossing effects can no longer be neglected, we can evaluate 
the scale down to which these reconstruction schemes can apply. 



can be understood partly from the fact that relatively few 
orders of perturbations theory are relevant at z — (since 
^s.c. = 9), so that it may not help much to explicitly include 
the effects of higher order terms. 

Reconstruction techniques are also used to sharpen the 
acoustic peak of the real-space correlation function or to 
restore the harmonics of the oscillations of the power spec- 
trum, in order to improve cosmological d istance measure- 
ment s and constrain ts on dark energy (|Eisenstein et al.l 
120071 : ISeo et al.l[2010h . Then, one can read in Table Q] the 
wavenumber k s . c . up to which one can hope to recover these 
baryon acoustic oscillations. At z = we can see that 
present schemes, which are based on the linear displace- 
ment field and manage to reach ~ 0.2h Mpc -1 (not neces- 
sarily for the amplitude but at least for the shape and po- 
sition of the oscillations) are not far from the upper bound, 
as could be expected from only a few orders of perturba- 
tion theory being relevant (n s c . — 9). At z = 3 it seems 
that one could go much beyond present schemes (which 
do not go much farther than 0.25/i Mpc -1 ), in agreement 
with n s c . = 66, which means that higher orders of pertur- 
bation theory are relevant. However, for the specific pur- 
pose of measuring cosmological distances from the baryon 
acoustic oscillations, the potential is limited by the relative 
amplitude of the oscillations of the linear power spectrum 
itself decreasing at higher fc, so that even a very good re- 
construction would not greatly enhance the signal-to-noise 
ratio. Nevertheless, pushing to higher orders (e.g., through 
resummation schemes) remains useful for other purposes, 
such as weak-lensing studies. 

3.3. Rise of power at the transition 

The increase of power on the transition scale to nonlinear- 
ity shown in Fig. [T] (especially in the left panel at z = 0) 
is reminiscent of a similar feature obser v ed for the gravita- 
tiona l dynamics (Hamilto n et al.l [l991l : iPeacock fc Doddsl 
Il996f ). This is usually interpreted from a Lagrangian point 
of view inspire d by th e spherical collapse dynamics. Thus, 
iPadmanabhanl (|1996h argues that, on these intermediate 
scales, one has £(x) oc £ L (q) 3 in real space with q 3 = 
x 3 (l+£). More generally, in Fourier space one writes for the 
nonlinear power per logarith mic interval of wavenum ber, 
A 2 (fc), the parametric system (jPeacock fc Doddslll996l ) 

k L = [l + A 2 (fc)r 1/3 fc, (47) 
A 2 (fc) = f[Al(k L )], (48) 

with a function / to be determined. These relations ex- 
press the conservation of matter, since the Lagrangian scale 
q ~ collapses down to the Eulerian scale x ~ 1/fc. The 
linear regime implies that f(x) ~ x for x <C 1, whereas 
in the highly n onlinear regime the stable-clustering ansatz 
([Peebles! [1982) gives the scaling f(x) - x 3/2 for x > 1. At 
the transition one obser yes a sharper growth, which is con- 
sistent with f(x) ~ a; 3 (jPadmanabhanl [l 996h . In practice, 
one builds a fitting formula for f(x) to match numerical 
simulations and to account for the dependence on the shape 
of the linear power spectrum. 

In any case, such models usually estimate the shape 
of the nonlinear two-point correlation function or of 
the nonlinear power spectrum by considering the col- 
lapse of a "typical" ove rdensity (jPadmanabhanl 119961: 
IValaeeas fc Schaefferl H997) (or merely obtaining f(x) from 
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Fig. 3. The nonlinear power per logarithmic interval of 
wavenumber, A 2 (fc), as a function of the linear power 
A|(fci), at the Lagrangian wavenumber defined by 
Ea. (j4"T)) . The lower curves are obtained from the Zeldovich 
power spectrum (|13|) , while the upper curves correspond to 
the power spectrum flT]) of our "sticky model" . The solid 
lines are for the redshift z = and the dashed lines for 
z = 2. 



simulations without further interpretation). It is interest- 
ing to note that this collapse also takes place within the 
Zeldovich dynamics studied here. In particular, the function 
6 = F{8l) that de scribes the sph erical collapse is no longer 
given by cycloids (|Peeblesl [l980') but by the simple expres- 
sion F(6l) = 1/(1 — ^l/3) 3 — 1 (|Bernardeau k, Kofmanl 
1995: IValageadl2009l) . (Collapse to a point is delayed from 
S c ~ 1.686 to S c = 3, since the motion does not accelerate 
as the gravitational potential well becomes deeper.) 

We show in Fig. [3] the functions A 2 (k) = f[A 2 L (k L )] 
defined by the system Il47|l - H8|) for the Zeldovich dynam- 
ics and the "sticky model". Indeed, from the knowledge 
of the nonlinear power A 2 (k) we obtain the Lagrangian 
wavenumber fcj, from Ea. (|47l) and next Aj(fci). For a given 
Eulerian wavenumber k, the Lagrangian wavenumbers kh 
obtained for both models, Zeldovich dynamics and "sticky 
model", are different. Then, Fig. [3] shows that, within the 
Zeldovich dynamics, this growth of the nonlinear density 
contrast through the spherical collapse is not sufficient to 
build up the increase in power on these mildly nonlinear 
scales, although one can see a modest rise at z = 2. Indeed, 
the shape of the nonlinear power spectrum is quickly gov- 
erned by the decay that takes place at higher k, due to the 
escape of particles beyond shell crossing that erases small- 
scale features, as recalled in Sect. 13.11 In contrast, such a 
growth is clearly seen in the "sticky model" . 

This means that such arguments, based on the evolution 
of a "typical" overdensity through the spherical collapse 
dynamics, are not good enough to explain the shape of the 
power spectrum on these scales since they clearly fail for 
the Zeldovich dynamics. Thus, one needs to consider the 
behavior of a whole range of typical density fluctuations 
(i.e., truly perform the Gaussian average over the initial 
conditions) , which includes a significant number of configu- 
rations where shell crossing has already taken place. Then, 
one cannot neglect the dependence on their behavior after 
shell crossing, and whereas the free escape associated with 



the Zeldovich dynamics is sufficient to erase any growth of 
power on the transition scale, a trapping of particles al- 
lows such a growth. In the "sticky model" studied here, 
this trapping after shell crossing is a simple sticking along 
one direction, as defined in (f3"4"|). which eventually leads to 
a universal k~ 2 tail at very high k for P s ticky(fc), as noticed 
in Sect. 13. ll In terms of the function /[A|(A;.l)], this corre- 
sponds to a large-a; scaling /sticky (x) oc 2; 3 /[ 2 ( n + 3 )] from the 
system igTfl - ggj) . 

In the gravitational case, this regime corresponds to the 
virialization within 3D bound structures, with an asymp- 
totic tail at high k that is poorly known. In particular, al- 
though overdensit ies seem to form ha los with an almost 
universal profile (jNavarro et al.l Il997l ). there is no first- 
principle derivation of the high-fc exponent of the den- 
sity power spectrum, and its degree of universality is not 
well known. However, especially in the case of CDM-likc 
power spectra that decay as Pl(/c) oc k~ 3 at high k, ha- 
los are rather smooth with a low amount of substructures. 
Then, the nonlinear logarithmic power spectrum does not 
seem to show a strong high-fc tail, such as the universal 
forms A 2 (fc) oc k 3 , oc k 2 or oc k 1 , associated with point- 
like masses, lines, or sheets, but appears to be consistent 
with a logarithmic or shallow slope asymptote ([Smith et al.l 
l2003t ). This means that the steep rise in the function 
A 2 (fc) = /[A 2 (&£,)] observed in the intermediate regime, 
1 < A 2 (fc) < 100, breaks down at A 2 (fc) ~ 100 to reach 
this slow growth at higher k, in contrast to the "sticky 
model" considered here where there is a regular transition 
from the linear regime, where /sticky (x) ~ x, to the highly 
nonlinear regime, where /sticky (x) oc x 3 ^ 2 ( n+3 ^ . 

On the other hand, in agreement with the fact that the 
linear density contrast associated with full spherical col- 
lapse to a point is delayed from 8 C — 1.686 to 3, the rise ob- 
tained in Fig. E] at A|(fc) ~ 1 appears at somewhat higher 
values of A 2 (k) than would be the case for the gravitational 
dynamics, as can be che cked by compariso n with the results 
of ./V-body simulations (|Smith et al.ll2003j) . In this respect, 
the spherical collapse captures some of the properties of the 
early rise of A 2 (fc). 

4. Redshift-space power spectrum 

The density field observed from galaxy surveys is distorted 
by their peculiar velocities, which i ntroduce an error in 
the m easure of their radial position (|Jacksonlll972t iKaiserl 
Il987f ). Indeed, both the mean Hubble flow and the peculiar 
radial velocity of a galaxy contribute to the redshift of the 
observed emission and absorption lines. On the other hand, 
the measure of the angular positions on the sky are not in- 
fluenced by these peculiar velocities. In a "plane-parallel" 
approximation, where the line of sight has a fixed direc- 
tion e z , radial and angular modes correspond to longitu- 
dinal and transverse wavenumbers k\\ and k^. Then, the 
power spectrum measured along the transverse direction 
(i.e. fc|| = 0) is not affected by these "redshift distortions" 
and is equal to the real-space power spectrum considered 
in the previous sections, whereas the power spectrum mea- 
sured along the longitudinal direction (i.e. kj_ = 0) is mod- 
ified. In particular, this means that the power spectrum is 
no longer isotropic. 

This implies that some additional information is con- 
tained in the redshift-space power spectrum, since by com- 
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paring the longitudinal and transverse components, we can 
derive some information on the velocity field and obtain fur- 
ther constraints on cosmology. For instance, within linear 
theory, comparing the redshift-s pace power spectrum over 
several directions (|Kaiserl fl987) or expanding its angular- 
dependence over Le gendre polynom ials and comparing the 
first few multipoles (jHamiltonl 1992) , one can constrain the 
ratio / of the linear growth rate to its value in a critical den- 
sity universe, hence ft m since / ~ ^m 6 - On the other hand, 
while clustering along the transverse directions of statisti- 
cal standard rulers such as BAO gives the angular distance, 
real-space clustering along the longitudinal direction gives 
the Hubble rate, so that complementary information can 
be derived from both directions, and one can use the lon- 
gitu dinal/transverse ratio to pe rform the Alcock-Paczynski 
test (jAlcock fc Paczvnskilll979] ). 

Thus, it is interesting to extend the analysis presented 
in the previous sections to the redshift-space power spec- 
trum. Rather than expanding on multipoles, we focus on 
the clustering along the radial direction, as compared with 
the transverse directions. Of course, another effect that 
comes into play in galaxy surveys is the bias, which may 
show some scale dependence. However, we do not study this 
effect here, as this is a rather different process. 



4.1. Nonlinear Zeldovich power spectrum in redshift space 

We first recall the nonlinear redshift-space power spectrum 
associated with the Zeldovich dynamics. The redshift-space 
coordinate s of a galaxy is 



s = x 



aH 



(49) 



where x is its comoving position, v = ox its peculiar veloc- 
ity, and e z the unit vector of the line of sight. Within the 
Zeldovich dynamics ((IJ the peculiar velocity is 



aD 
IT 



leading to 

s(q, t) = q + * L + / (e z ■ * L ) e z 



(50) 



(51) 



where D is the linear growth factor and f(z) = 
dlnP/dlna. In the following we use a "plane-parallel" ap- 
proximation and we denote the longitudinal and transverse 
directions to the line of sight by the subscripts || and _L 
Thus we have 



3|| =g|| + (l + /)*£]), S ± =q ± + * L _L. 



(52) 



The conservation of matter reads again as p s (s)ds = pdq, 
where we denote the redshift-space quantities by a super- 
script "s", and as in Eqs.Q-® the redshift-space power 
spectrum reads as 



P s (k) 



dq 

(2tt) s 



o ik-[s(q)-s(0)h 



Following iTavlor fc Hamilton! ()1996D . it is convenient to in- 
troduce the vector K, which is the wavevector k stretched 
by (1 + /) along the line of sight, 



K=(l + /)fc|, e z +k ± . 



(54) 



Then Eq.(|55]l writes as (see Ea.(fT3"|l) 

ps(k) = I dq c ik-q.( c lK-r*z,(q)-*tro^ 



(2tt) s 
(2tt) 3 



kq -if 



V?-io(9)-(l-3/4 ? )i2(«)] ; (55) 



where /iKq = (K • q)/ (K q) . Using spherical coordinates 
about the vector K, expanding parts of the exponentials 
and using Eo. (fl4"f we obtain 



P s (k) 



dqq 2 
„ ~2^~ 



e -K 2 [al-I (q)+2I 2 (q)] 

( 

{6K%(q)\ e (2kq(nl k -l) 



£,m=0 



{l + m)\ 
l\ [(2m) 



1112 



V kq/j, 



with 



k • e. 



i+k = 



Kk 



HKk = 



\ ^Kk 



jt+m{kqiJLKk), (56) 



K k 



(57) 



K 2 = k 2 [l + ^ k (2f + f)}. (58) 

Expression (|56p holds for any wavenumber k and shows 
how the redshift-space power spectrum depends on both k 
and /j,k- As recalled above and as is obvious from Eq. (|52"l) . 
the redshift-space power spectrum for wavevectors perpen- 
dicular to the line of sight is equal to the real-space power 
spectrum, 



Pl(fe ± ) = P s (k ± ) =P(k±) 



(59) 



and we can check that we recover Eq. (fT5")) from Eq. ([5rj|) 
for /Xfc = 0. For longitudinal wavevectors, we simply have 
K = (1 + /)k, so that Eg. ([55]) simplifies as 

fork^ie,: Rf (*„) = P s (k\\e z ) 



dC > „ik q „-fc 2 (l+/) 2 [^-/o(9)-(l-3p 2 )/2(9)] 

(2tt)3 



(60) 



where \x = (k • q)/(fcq) = q\\/q, as in Eg. (fl~5)) . Thus we see 
at once that, along the line of sight, the effect of redshift 
distortions is merely to multiply the amplitude of the linear 
power spectrum by a factor (1 + /) 2 ; that is, we can absorb 
the factor (1 + f) 2 into o 2 , Iq and I2, or more simply into 
Pi(fc), as was obvious from Eg. (152")) . Therefore, we directly 
obtain from Eg. (|15[) the expression 



Jo 



dqq 2 , 



2tt 2 



k 2 (l+f) 2 [ai-I (q)+2I 2 (q)] 



(61) 



since we only need to multiply each term that involves 
Pl(/c) by a factor (1 + /) 2 . In a similar fashion, the standard 
perturbative expansion (|19l) becomes 



(53) Pf(fc) = J2 P \\ ( fc ) with P || ( fc ) = (l+f) 2n P {n) (k), (62) 



whereas the "renormalized" perturbative expansion (|24p 
becomes 



P|(fc) = e- fe ^ 1 +^^p£ ) (fc), 



(63) 



n=l 
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with p£\k) = (l + f) 2n P£\k). 



(64) 4.3. Redshift-space logarithmic power 



Thus, each order n of the perturbative expansions gets mul- 
tiplied by a factor (1 + f) 2n as we go from real space to 
redshift space, for modes that are aligned with the line 
of sight. In particular, at lowest order we have -Pif(fe) = 

(1 + f) 2 PL(k) + so t hat we recover the well-known boost 
factor of Kaiser (1987), associated with the infall of galax- 
ies within gravitational potential wells. On the other hand, 
the factor (1 + /) 2 also sharpens the Gaussian damping 
prefactor in the expansion (|63p. 



4.2. "Sticky model" nonperturbative correction 

We now turn to the nonperturbative correction associated 
with the "sticky model" introduced in Sect. 12.41 Again 
we focus on the redshift-space power spectrum Pif(fc) for 
wavenumbers that are parallel to the line of sight. From 
Ea. ([55|) we now need to compute the average (e lk As ) q . 
Before shell crossing, or more precisely, as long as Axli > 
for a Lagrangian-space separation vector along the ei-axis, 
the "sticky model" follows the Zeldovich dynamics, 

k = ke z , q = qe x : 

if A* L1 > —q : k • As = k ■ q + (1 + /)k ■ * L , (65) 

whereas once Ax^i < we affect the same position, Axi — 
0, and the same peculiar velocity, Av-y = 0, along the ei- 
axis, for both particles as in (|34p. whence 



if A* L i < —q : As\ = 0, hence 

k.A B =(l + /)[fc2*£2 + fc3*L3]. (66) 

(The line-of-sight axis e z and the axis ei of the Lagrangian 
separation vector q are not related.) At the shell-crossing 
time we have a discontinuity in the redshift-space separa- 
tion As, since Asi jumps from —fq to zero. Indeed, just 
before contact, the two particles move closer, with the finite 
relative velocity Avi — —(aD/D)q, while after collision 
their relative velocity is set to zero. This should actually 
be understood in a loose sense, since within the "sticky 
model" we consider neither transverse directions nor the 
cloud-in-cloud problem. Then, as in Sect. 12.41 the redshift- 
space power spectrum of the "sticky model" reads as 



P/ ticky (k)=P| el (k) + ^c.(k), 



(67) 



and from Eqs. (|65p ~ (|66[) we obtain for wavenumbers along 
the line of sight 



dqq 2 _„2/,n^ 



p> (*)= / ^e-' / dpe-w+n t 1 -™*- 10 - 1 '! 

e -*fw w [ \_ J JH II I . ( 68 ) 



x Re < w\ 



y/2a\\ 



V2a« 



As compared with Ed.(|45[). the new exponential factor 
g-ifc/w j s (j ue ^ ^fie discontinuity of As at AvP^i = —q, as- 
sociated with shell crossing within the "sticky model" . The 
factors (1 + /) that multiply the longitudinal wavenumber 
k could be expected from Eas.(|65 p -(|66 p . 



As for the real-space Fig.[TJ we show in Fig.H]our numerical 
results for the redshift-space logarithmic power, for longi- 
tudinal wavenumbers k, defined as 



Af(k) = 47rfc 3 P | f(/c). 



(69) 



We use the same definition (j4"6"]l . even though Pu(k) only 

holds along the longitudinal direction and PJ(fe)(= P(k)) 
holds along the two transverse directions (so that using a 
factor k instead of k 3 would be more natural here), to make 
the comparison with Fig. Q] easier. In particular, this means 
that Aj s (fc) is given by Fig. [T] for transverse wavenumbers. 

As compared with Fig. [TJ we can check that the linear- 
regime power, f || j is amplified by a factor (1 + /) 2 , 
whereas standard higher order terms are amplified by fac- 
tors (1 + f) 2n , see Eq. (j62]) . Even though the same factors 
apply to the higher order terms of the "renormalized" per- 
turbative expansion, see Eq. (|64l) . their peak height is not 
greater than in Fig. [T] for n > 2 because of the stronger 
Gaussian damping prefactor e~ k "v in Eg. (1551) . For 

large n this even leads to a smaller amplitude as compared 
with the associated real-space contribution. This makes the 
full nonlinear Zeldovich and "sticky model" power spectra 
greater than their real-space counterparts on large scales, 
in the quasi-linear regime, but smaller in the highly nonlin- 
ear regime. In particular, it is clear from Eg. ([50)1 that the 
high-fc damping of the Zeldovich power spectrum becomes 
sharper because of the factor (1 + f) 2 in the exponential 
term (even though this only leads to a power-law decay as 
noticed in Sect. 13.11) . As well as for the real-space power 
spectrum, the "renormalized" perturbative expansion (f6"3"f 
is more convenient to distinguish the relative contributions 
of higher order terms, as they are all positive and do not 
show the cancellations associated with the standard expan- 
sion ([52]). 

Since within the "sticky model" we have after shell 
crossing Asi = 0, just as we had Ax\ = in real space, 
we again obtain a k~ 2 tail at high k for the nonperturba- 



tive term P, 



(k), whence A| s s t . cky (fc) 



k, as would be 



the case for planar structures in redshift space. However, 
because the dominant effect of redshift distortions on small 
scales is to decrease the power, through damping factors of 
the form e ( L +J) a v j the nonperturbative redshift-space 
contribution P| s c (k) is less than its real-space counterpart. 

Coupled with the faster decay of P| Zc i(fc), this leads to a 
temporary decrease of A 2 ^ ticky (k) at k ~ 1 to Ah Mpc -1 at 
z = 2, before the asymptotic tail oc k becomes dominant. 

The "sticky model" is not intended here to describe the 
redshift-space power spectrum better than the Zeldovich 
approximation. Indeed, on small scales setting Avi = 
is not realistic, since one should rather describe multi- 
streaming and take into account the finite velocity disper- 
sion of collapsed objects built by virialization processes. 
Then, instead of a "sticking" in redshift space one obtains 
a finite size extension that is greater than the real-space 
size of the object (because of the additional term associ- 
ated with the velocity dispersion in Ea. (|4"9)l ). As is known, 
this explains th e characteristic "fingers of god" ob served in 
galaxy surveys (|Jacksonlll972t iTadros et al.lll999l) . and in- 
stead of a universal growing tail A|g ticky (fc) cx k, this leads 
to a decaying tail (as for the Zeldovich power spectrum but 
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Fig. 4. The logarithmic redshift-space power spectrum, A| s (fc) = 47rfc 3 P|j i (fc), for wavenumbers along the line of sight. 

Both plots, at z = (left panel) and z — 2 (right panel), are the redshift-space counterparts of those shown in Fig. [T] for 
the real-space power spectrum. 



with different quantitative properties, since the latter are 
set by the complex nonlinear virialization processes rather 
than the linear velocity field). Of course, it is possible to 
modify the "sticky model" to implement such a velocity 
dispersion and to obtain more realistic results in redshift 
space. For instance, while retaining "sticking" in real space, 
one could allocate a random relative velocity to the pair, 
so that the perturbative expansions remain unchanged and 
the nonperturbative contribution involves a damping term 
of the form e~ fc where ovir is the characteristic velocity 
dispersion in virialized objects. However, this goes beyond 
the scope of this article, since here we are only interested in 
the comparison of perturbative and nonperturbative terms 
and the simple "sticky model" in the form (|66l) is sufficient 
to estimate the scale where nonperturbative corrections be- 
come important. 

4.4. Redshift-space perturbative and nonperturbative 
contributions 

Finally, as for the real-space Fig. [51 we show in Fig. [5] the 
ratios of the successive terms of the "renormalized" pertur- 
bative expansion (|63p with respect to the nonlinear power 
spectrum of the "sticky model" , for the same orders n = 1 
to 5, as well as n = 10, 15, 20, 30, 50, and 70. We also plot 
the ratio P| s c (fc)/-P|| sticky (^)i to compare with the am- 
plitude of the nonperturbative correction associated with 
shell crossing effects, and we consider the four redshifts 
z = 0,1,2, and 3. As for the real-space power spectrum, 
the potential of perturbative expansions is greater at higher 
redshift, as we can include more perturbative terms before 
they become subdominant as compared with the nonpertur- 
bative contribution P,? c (k) , and they extend by a larger 
factor the range of k where systematic results can be ob- 
tained. 

We give in Table [5] the wavenumbers fcy lloop and feij s c 
where the second order of perturbation theory, P^\k), 
and the nonperturbative contribution, Pf (fc), reach levels 



z 


% 


fcffiioop [/iMpc- 1 ] 


fejf.. c . [ftMpc" 1 ] 






1% 


0.025 


0.13 







10% 


0.052 


0.26 


8 




50% 




0.57 






1% 


0.03 


0.3 




1 


10% 


0.064 


0.63 


18 




50% 




1.3 






1% 


0.033 


0.68 




2 


10% 


0.08 


1.5 


37 




50% 




4 






1% 


0.043 


1.3 




3 


10% 


0.1 


3.2 


72 




50% 




9.8 





Table 2. The dependence on redshift z of the wavenum- 
bers where the ratios e~ k2{1+f)2 ^ (k)/P^ sticky (k) (at 

*ifuoop) and ^Is.c.W/^fstickyW fcf s .J reach 1%, 10%, 
and 50%. The last column gives the last order njj s c of the 
"renormalized" perturbative expansion that is not fully be- 
low the nonperturbative term. This table is the redshift- 
space counterpart of the real-space Table [TJ for longitudinal 
wavenumbers. 



of 1%, 10%, or 50%. As noticed above from the comparison 
of Figs. [5] and [21 the interval [&jj lloop , fejj s c ], where higher or- 
ders of perturbation theory apply, shifts to higher k and be- 
comes broader at higher redshift, as for the real-space spec- 
trum described in Table [TJ However, it is shifted towards 
smaller k than with the real space power spectrum, as could 
be expected from the fact that redshift-space perturbative 
terms of order n are multiplied by prefactors (1 + f) 2n , 
see Eas. ([(32"j) - (|f)4")) . so that higher order terms become rel- 
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Fig. 5. The ratios of the "renormalized" perturbative terms of the expansion (l63l) to the nonlinear "sticky model" 
redshift-space power spectrum, and of the nonperturbative contribution Pfi sc (k). This is the redshift-space counterpart 
of the real-space Fig. [SJ for longitudinal wavenumbers. 



evant earlier and on larger scales. On the other hand, the 
order nf, a c , after which perturbative contributions become 
subdominant as compared with the nonperturbative con- 
tribution, is similar to the one obtained for the real-space 
power spectrum, and grows from 8 at z — up to 72 at 
z = 3. 

Since redshift distortions are greatest along the line 
of sight and they vanish for transverse wavenumbers, one 
would obtain similar results for multipoles of the redshift- 
space power spectrum, where one expands the dependence 
on the angle with respect to the line of sight of P s (k) over 
Legendre polynomials of \x = (k ■ e z )/k. On a quantitative 
level, one would obtain results in-between those presented 
in Tables Q]and H 

5. Conclusion 

In this article we have introduced a "sticky model" that 
coincides with the usual Zeldovich dynamics before shell 
crossing, while after shell crossing it includes a sticking of 
particle pairs along their longitudinal direction. This im- 
plies that the nonlinear density power spectra of both mod- 
els have the same perturbative expansions and only differ 
by nonperturbative terms that arise from the dynamics be- 
yond shell crossing. Since we can obtain explicit expressions 
for perturbative terms at all orders and for this nonper- 



turbative correction, we have been able to compare their 
respective amplitudes in detail, focusing on a ACDM cos- 
mology. 

In real space, we find that perturbative approaches 
based on the fluid description can extend the wavenum- 
ber up to which systematic analytical predictions can be 
obtained from k ~ 0.033ft Mpc -1 up to ~ 0.23/i Mpc" 1 
at z = 0, and from ~ 0.07ft Mpc" 1 up to ~ 2.2ft Mpc" 1 
at z = 3, as compared with linear theory (if we need an 
accuracy of 1%). We also give detailed results for other 
redshifts and for 10% and 50% accuracy levels. Going to 
higher k requires taking shell-crossing effects into account. 
Nevertheless, these results show that such perturbative ap- 
proaches, based on hydrodynamical equations of motion, 
have a significant potential, especially at high and moder- 
ate redshifts. Since it is unlikely that it will be possible to 
explicitly compute high-order perturbative terms up to the 
last order n s c . which is above the nonperturbative correc- 
tion (n s . c . = 9 at z — and n s . c , = 66 at z — 3), this pro- 
vides a strong incentive to develop resummation schemes. 
Then, the hope is that such methods can efficiently re- 
sum most of the contributions associated with higher order 
terms and achieve a faster convergence. 

This analysis is also useful in the context of cosmologi- 
cal reconstruction, where one tries to recover the dynamical 
history of a given region of the sky from its present den- 
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sity field (and next to derive from this reconstruction the 
present velocity field). Indeed, these methods usually ne- 
glect shell-crossing effects. In particular, we find that the 
Monge-Ampere-Kantorovich method is close to optimal at 
z = 0, as it fares well down to the scale where shell-crossing 
contributions dominate. With respect to the reconstruction 
of the baryon acoustic peak or of the acoustic oscillations of 
the power spectrum, there seems to be only modest room 
for improvement over current methods based on the linear 
displacement field at z — 0, in agreement with only a few 
orders of perturbation theory being relevant (n s .c. = 9). At 
higher z one could in principle do much better by including 
higher order terms, but since baryon acoustic oscillations of 
the linear power spectrum have a small amplitude at high fc, 
this may not significantly improve the signal-to-noise ratio. 
However, this remains useful for other observational probes 
that are sensitive to the shape and amplitude of the power 
spectrum on weakly nonlinear scales, such as weak-lensing 
surveys. 

We have also pointed out that the behavior of the sys- 
tem after shell crossing plays a key role in the shape of 
the density power spectrum on mildly nonlinear scales, 
1 < A 2 (fc) < 100. Thus, arguments based on the spherical 
collapse dynamics are not sufficient to explain the steep rise 
of the nonlinear power in this intermediate regime. Indeed, 
while for the Zeldovich dynamics such a growth is almost 
entirely wiped out by the escape of particles to infinity (as 
the random linear displacements erase small-scale features 
and lead to a decay of A 2 (fc) at high k), for the "sticky 
model" there is a regular growth up to the high-fc asymp- 
tote, P s ticky(^) ~ k~ 2 . Therefore, it is not sufficient to con- 
sider the spherical collapse of a "typical" overdensity to ex- 
plain the shape of P(k) on these scales, as one must take the 
Gaussian average over a whole range of initial density fluc- 
tuations, which show a significant amount of shell crossing 
that has a strong impact on the resulting power spectrum. 

Finally, we have obtained similar results for the redshift- 
space power spectrum. Again, the scope of perturbative ap- 
proaches (based on the fluid description) is greater at higher 
z. As compared with the real-space power spectrum, the 
characteristic wavenumbers where higher order terms of the 
perturbative expansions and nonperturbative contributions 
come into play are shifted to lower values of k. Indeed, be- 
cause of the amplification of perturbations from the uniform 
Hubble background by the additional contribution from pe- 
culiar velocities, the nonlinear regime is reached on larger 
scales than in real space. 

The "sticky model" presented in this article, and more 
generally the expression ([7]) of the density power spectrum 
in terms of the Lagrangian displacement, could serve as a 
basis for models for the nonlinear power spectrum. In par- 
ticular, the spirit of this model could be incorporated into 
the Lagrangian-based models, associated with the paramet- 
ric system (jT7)) - (P5| . or into the "halo model", to make 
the bridge between the perturbativ e regime and the highly 
nonli near, nonperturbative regime (Valagea s fc Nishimichil 
l2010f) . This could also be extended to redshift space by in- 
cluding a model for the velocity dispersion within virialized 
objects. However, this is beyond the scope of the present 
article so we leave it to future works. 
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